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ABSTRACT 

Analytical expressions are presented for relativistic corrections to the Kompaneets 

equation and the Sunyaev-Zel'dovich y-distortion. The latter provides a convenient 

method of inferring both the temperature and Thomson optical depth of the gas in 

clusters of galaxies using only observations of the spectrum of CMBR photons passing 

through the cluster. The relativistic correction gives an additional component of the 

Q\ ■ S-Z effect apart from the y-distortion and kinematic S-Z effect. The perturbative 
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Q^ ■ but not for higher temperatures. 
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method used is shown to provide a very accurate approximations when kT e < 100 keV 



Subject headings: radiative transfer — scattering — cosmic microwave background 
> ■ galaxies: clusters: general 



1. Introduction 

In the seminal paper, Zel'dovich and Sunyaev (1969), it is shown how interactions of 
primordial photons with a hot plasma will lead to deviations from a blackbody spectrum of 
photons. The classical Sunyaev-Zel'dovich (S-Z) effect gives the first order effect on the spectrum 
as one increases the temperature of the electrons from zero and is not valid as the electron 
velocities become relativistic. In this Letter we give the next order effect, which is important 
for hot gas in clusters of galaxies, describe a perturbative expansion for computing higher order 
corrections, and consider the convergence of this expansion. 



The classical S-Z distortion was based on the Kompaneets equation ( Kompaneets 1957| ) 



which describes how the spectrum of photons evolves under the action of Compton scattering off 
of a stationary plasma of hot electrons. To a first approximation, for non-relativistic electrons 
(v 2 <C c 2 ) and soft photons (e <C m e c 2 ), kinematics dictates that Compton scattering does not 
change the photon energy and there should be no change in the photon energy distribution. The 

2 

Kompaneets equation describes the rate of change in the photon spectrum to order ^ and ^ . 
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The next order term, 



oc 



or more precisely oc 



\m c c J 



, was presented in |Stebbins 1997a| and 



in this Letter we summarize these results while adding other analytical results and applications. 
Compton scattering has been well studied, by a variety of techniques and non-perturbative 
computations of the spectral distortion of an incident blackbody spectrum passing through a hot 
gas at arbitrary temperatures have been published (Fabbri 1981, [Rcphaeli 1995 ), The perturbative 
results of this Letter complement the more precise numerical results in that they are analytical 
and describe the corrections to the classical y-distortion in terms of a given spectral shape. As we 
shall see, this next order correction provides an good approximation to the numerical results in 
the temperature regime relevant to the hot gas in clusters of galaxies. Rephaeli and Yankovitch 
(1997) have already shown that these corrections, combined with relativistic corrections to the 
X-ray emission can have significant effects on the estimates of the Hubble constant from cluster 
observations. A more expansive discussion of the perturbative results presented here may be found 
in Stebbins (1997a,b). 



2. Collisional Boltzmann Equation 



We describe the state of the primeval gas of photons in terms of the quantum mechanical 
occupation number in phase space, n 7 . The evolution of n 7 can be described by the collisional 
Boltzmann equation which has the form 



Z)n 7 (p 7 ) 



Dt 



C( Po 



(1) 



where C(p 7 ) is the scattering term which describes the effect of scattering and A is a convective 
derivative along the photon's trajectory in phase space. Assuming isotropy of the photon and 
electron distribution function, we need only solve for the change in n 7 as a function of the photon 
energy, e, and not momentum, p 7 .[] The collision term for Compton scattering may be written 
(Stebbins 1997a,b) 



C(e,A) = J d(3f((3) J dA 



' ;S( v ^,A)(l + n~ f (e))n^ 



-5(e,A)(l + n 7 (e(l + A)))n 7 (e) 



(2) 



where S(e, A) gives an angle-average of the differential cross-section: 

S(e,A) = cN c ((l-fa)^ . 



(3) 



If the photon are not at rest with respect to the electrons then their distribution functions will not both be 
isotropic in the same rest frame, which will lead to a direciion-dependent spectral distortion known as the kinematic 
Sunyaev-Zel'dovich effect, which is a type of Doppler shift. 
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Here N e is the electron density, n is the incident electron direction, (3 is the electron velocity in 
units of c, f(8) gives the distribution of electron speeds, e is the incident photon energy, and A is 
gives the fractional change in the scattered photon, i.e. e' = e(l + A). The 1st term in the collision 
integral gives the scattering into the beam and the 2nd term the scattering out of the beam. 



3. Fokker-Planck Expansion 

As previously mentioned, for non-relativistic electrons and soft-photons, kinematics dictates 
that the change in the photon energy is small so that the dependence of 5(e, A) on it's second 
argument will be very sharply peaked around A = 0. The rest of the dependence on A in eq. [2] 
is much smoother, and may be approximated by it's Taylor series about A = 0. This Taylor 
expansion yield a Fokker-Planck type of approximation. Taylor expanding the A-dependence of 
the distribution function n 7 leaves only derivatives of n 7 evaluated at e. These derivatives may 
be taken out of the A and (5 integrals and one is left with a differential equation for n 7 rather 
than an integro-differential equation. The coefficients in the differential equation are determined 
by the moments of the A distribution, which can be determined analytically when the electron 
distribution is thermal. 



The full Klein-Nishina cross-section is a rather complicated function, and for the purposes of 
the CMBR we are mostly interested in the soft-photon limit. This leads one to Taylor expand the 
full cross-section in powers of a = and at each order compute eq [| The 1st term is (Stebbins 
1997a,b) 

S(e,A) = cN c a r f d/3 3F(A, /3sgn(A)) + 0{a l ) (4) 
Jo 



where 



F(A, b) = sgn(A) x H(l - (1 ^ ) x 



3(1 - b 2 ) 2 (3 - b 2 )(2 + A) (1 - 6)(1 + A) 



166 6 



l + b 



+ 3(1 ^(1+% b)A) (4(3 " + + 2(6 + b " ~ b3 + 2tA)A + (1 " &2)(1 + b)A2] 



(5) 



where TL is the Lorentz-Heaviside function. With explicit functions such as these it is tedious but 
straightforward to compute the Fokker-Planck (=Taylor series) expansion of eq ^ in powers of a 
and B or for a thermal electron distribution in powers of O e = % . Much of the methodology for 
such expansions was developed in Barbosa (1982). 

One may write this Fokker-Planck expansion as 



~d7 



]_d_ 

e 2 de 



(n,m) 
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(6) 
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where /T (n ' m) [n 7 ] indicates a, potentially nonlinear, differential operator acting on n 7 which does 
not depend on m e c 2 . The pre-operator, 4rjj, is guaranteed by the form of eq |2| which explicitly 
conserves photon number. We have used the Thomson optical depth: r = / dt N e aT c where <tt is 
the Thomson cross-section. One finds that to include all terms that contribute to K {n ' m) one must 
Taylor expand to order 2n + m in A. Since the energy redistribution is a relativistic effect we find 
that K^°' 0) = 0. The classical Kompaneets equation contains the two terms: 

K W) [ nj ] = e K W [ n7 ] = (i + n7 ) ^ . ( 7 ) 

The 2nd order A-expansion of Barbosa (1982) is only accurate enough to infer the one additional 
term 

K<°' 2 >[n 7 ] = - — (l + n 7 )n 7 . (8) 

The new result (Stebbins97a,b) is 

,,,, m 5 dn-y 21 9 d 2 n^ 7 q 9 3 n^ . . 

which provides the next order correction for hotter electrons, but not for more energetic photons. 
Challinor & Lasenby (1997) have confirmed these corrections using an independent derivation and 
have also computed some additional terms. Combining this new terms with the classical result we 
find an extended form of the Kompaneets equation: 



<9n 7 1 d 
dr e 2 de 



., / kT c / 5_ \ dn-y 7 ^ / 9 d 2 n~ t , d 3 n'. 
e % 1 + 7^@c e ^r 2 + T7T©o 6e 2 — ^ + e 3 — ^ + a 1 + n 7 n n 



m P c* 



2 1T 2" e 5e ^1() e T cV ^ c 5e 3 -r«v--r-r/»7 



(10) 



which will be more accurate at higher electron temperatures than the classical form. 



4. Extended Sunyaev-Zel'dovich Distortion 

The idea of the S-Z distortion is that one starts out with a background radiation which is 
close to a blackbody spectrum, just what we expect to be produced by the early universe, and it is 
slightly distorted by the action of hot ionized gas through the Compton scattering process we have 
just described. In this small distortion limit we need just substitute a blackbody spectrum with 
temperature T 7 for n 7 in the right-hand-side of the Kompaneets equation and integrate over r to 
obtain the distortion to the spectrum. If we apply this technique to the Fokker-Planck expansion 
we obtain a series for the total distortion: 

/kT 
drQ^Q™ 6 7 = — \ (11) 
ii < i in 1 1 TTl e C 

where 

A^r' (*) = K <-») [nf] < B = ^— T x = ^ . (12) 
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From the known K (n ' m) one finds 
Anff(x) = 

An™(x) = -A„"(x) = ±A*™(x) = (x J±l - 4) 

. , am/ . xe x / 47 e x + l 42 2 e 2:c + 4e a; + l 7 = (e x + l)(e 2:c + lOe* + 1 
An 2 7 0) x = ^ -10 + —x x 2 — -75— + —x 3 v 



(e 8 - l) 2 \ 2 e x - 1 5 (e 1 - l) 2 10 {e x - l) 3 J 

(13) 

Since one expects that to each order in energy that a blackbody spectrum is a stationary solution 
when the electron and photon temperature are equal, one expects the sum rule 

N 

^An^(x) = 0. (14) 

n=0 

This sum rule tell us that between the N + 1 terms in this sum only N give linearly independent 
spectral shapes for the distortion. Applying the sum rule to iV = 2 one can infer the function 
Angz^(:r) via the relation Ang 1 ^ = -An^ 01 — An s °z 2) (x). 

The classical S-Z y-distortion includes only An S z 0) and An!^ 1 ' and is parameterized by the 
"y-parameter" : 

y = y(M) _ yW.) = f dT k ^ ~ T 7) (15) 
J m e c z 

Additional terms in the Fokker-Planck expansion introduces additional parameters such as Y^ 2 ,o) 
which are measures of different moments of the electron or photon temperature. 



5. Accuracy and Convergence 

Expansions like those of eq.s Pfe|i~T| are most useful when good accuracy is obtained with only 
the first few terms. These expansions give the correct Taylor series expansion in e and 7 about 
zero, but this may not yield a good approximation even when these variables are small, and in fact 
these series may not converge for any value of the parameters. As a first check we have compared 
the series expansion of An 7 to order and O 2 to the numerical results of Rephaeli (1995) in 
fig H for a = and kT c = 5keV, lOkeV , 15keV (9 C = 0.0098, 0.0196, and 0.0294). The truncated 
series yields a good approximation to the numerical results, but with discernible differences. While 
it is perhaps surprising that such large corrections to the classical S-Z distortion are present at 
such low temperatures, the fact that the next order correction makes up for most of the total 
correction indicates that this series will be useful. 

To examine the accuracy and convergence of the series expansions one may look at the low 
frequency limit of the S-Z distortion. By substituting the power series expansion in x of the 
blackbody spectrum into eq[2] one obtains a power series expansion for the extended S-Z distortion: 

Anf(x) =r_ 1 (0 c )i+r 1 (e c )x + r 3 (e c )x 3 + 0(x 5 ,ei) (16) 
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Fig. 1. — Plotted vs. frequency is the shape of the spectral distortion (x 3 An(x) oc brightness) 
to thermal CMBR photons passing through hot gas. Shown are three triplets of curves for gas 
temperatures of 5, 10, and 15keV in order of increasing amplitude. In each triplet the dotted curve 
shows the classical y-distortion, the dashed curve the 0(0;?) perturbative prediction, and the solid 
curve the "true" distortion taken from the numerical results of Rephaeli (1995). 
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the even powers having zero coefficients. Here we will be interested in the "Rayleigh-Jeans" term, 
r_i(O e ), which gives the fractional change in brightness at the lowest frequencies. This is given by 



r_i(G c ) = / dp- 



o 



9 - 15 13 2 + 10/? 4 3(3-6/? 2 + 4/? 4 -/? 6 ) 1 1 + /3 

m 



4/3 4 4/3 5 1-/3 

+ 0(0*) (17) 



3(3-/? 2 )(l-/? 2 ) 2 f h } + ^' 2 



16/3 6 V 1-/3 

where 7 = (1 — /3 2 )~2 and -K2 is a modified Bessel function. While one may not be able to express 
this integral in terms of simple functions it is easily evaluated numerically. One may also compute 
the Taylor series of r_i(0 e ) about zero: 

„ l( e o) = - 2Sc + ^ - ™e> + 32»ej - + + o(@l x \ e>). (is) 

11 ej 5 20 140 320 112 c e 7 



which also gives the low frequency limit of the expansion of eq 11 with m = 0. The convergence 
of this series should tell us something about the convergence of eq 11. In fig || is plotted the 
fractional error of truncated Taylor series approximations to r_i. It seems clear that the series 
does not converge for kT e > 80keV and likely that it has zero radius of convergence. However 
below ~ 100 keV the 1st few terms in the Taylor series do provide excellent approximations to r_i. 
If we take this to be representative of the more general series of eq |ll| we may conclude that these 
extensions to the S-Z distortion should be sufficient for application to clusters of galaxies where 
we expect kT e < 20keV. One just need be careful not to expand to such larger orders that the 
series starts to diverge.0 

It is not clear whether one might be able to slightly modify eq |ll] to obtain a series which 
is useful at all temperatures. We know that for ultra-relativistic electrons that essentially every 
scattering will take a photon out of the Rayleigh-Jeans region to much higher energies. Thus we 
expect the brightness at low frequencies to be degraded to a gray-body spectrum, the brightness 
decreased by the factor e~ T . This reasoning correctly predicts the limit r_i(+oo) = —1. This 
asymptotic behaviour makes the error of any approximant which is a polynomial in e diverge 
at large e - Once can hope to improve convergence if one instead Taylor expands in some other 
variable X = f(Q e ) where /(0) = and /(oo) is finite. For example if one takes X = O c /(l + 30 e ) 
the Taylor series of r_ 1 in X, 

v ; 54 35 2240 448 2560 v ' ' v ' 

seems to have good convergence for all temperatures. Note that the Taylor series in X contains 
no more information than that in O e , since one can derive one from the other. What is unclear is 
whether one can make something like this work for all frequencies. 



2 These convergence properties are indicative of an asymptotic series, an example of which is perturbation theory 
in quantum field theory. 
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Fig. 2. — Plotted versus electron temperature is the fractional error made by truncated Taylor series 
approximations to the Rayleigh- Jeans spectral distortion. The curves are: (dotted) the classical S-Z 
effect - a 0(0*) truncation; (short- dashed) a 0(Og) truncation; (long-dashed) a 0(@^) truncation , 
and so on up to an O(0 10 ) truncation. The hashed region to the right of the vertical line, and some 
of the hashed region to its left are accuracies which are not attainable at any order in the Taylor 
series. Some of the hashed region to the left of the vertical line is attainable by truncating at orders 
greater than C(6g°). Note that while the Taylor series is never very accurate when kT c > 100 keV, 
a truncated Taylor series can be extremely accurate at lower temperatures. 
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At cluster temperatures the accuracy obtained from the 1st few terms in eq [Tl| is so good 
that other effects may be the leading source of errors. One is the finite optical depth effect. If one 
expands in r as well as in temperature the leading r 2 term is that given in Fabbri (1981) which is 
~ 0(t 2 O 2 ). For clusters, where r ~ 0.01, this is liable to be much smaller than the y (2:0 ) An? 2 z 0) 
correction given above but should be roughly the same order as the next term, Y" (3i0 ) Are^ 3 z 0) , which 
we have not computed. To properly compute finite optical depth corrections in a non-uniform 
medium (such as a cluster) one must really do a proper radiative transfer computation. 



6. Application to Clusters 

Given that only the first few terms in eq |ll] give an excellent approximation to the exact result 
for cluster temperatures, we may think of the extended S-Z distortion of the CMBR spectrum 
as consisting a few distinct types of distortions corresponding to the different n and m in that 
equation. For the CMBR passing through clusters the terms with m ^ are negligible and we 
henceforth ignore them. Of the remaining terms the largest is An? z 0) (the classical S-Z effect) the 
2nd largest is An^ 2 z 0) , etc. The amplitude of the classical S-Z distortion is given by y = Y {1S)) and 
tells one about the product of the temperature and the optical depth. The amplitude of the next 
term, Y^ 2 ,o), can be combined with Y {lfi) to give a weighted measure of the temperature and optical 
depth separately 

kT-rnc^-l^^l ^_Y^ 2 _ Udr{kT c )f 

e cC y^°> " JdrkT e ~ Jdr(kT c ) 2 ' { U) 

If the electron temperature was constant along a line-of-sight, then we expect each of the ratios, 
^Wn.o') ^ ) to be the same independent of n. However one can use deviations from this rule to detect 

non-isothermality. For g clS ctS hot as 15keV ^£1 ~ 0.03 but since in the 100-400 GHz region 

An; 32 

gz ~ 5 the signal from the next order correction to the classical y-distortion is only about 7 

(1,0) 

times smaller than the classical effect. For many hot clusters one can expect this next order effect 
to be at the ^r- ~ 10~ 5 level which is certainly accessible with existing telescopes and should also 
be measurable by the all-sky survey of the PLANCK satellite. 

In order to measure quantities like Y^,o) one must be able to distinguish the An^ z 0) distortion 
from other effects such as anisotropics, the classical y-distortion, and contamination by dust, 
synchrotron, and free-free emission. In fig [|| we plot the shape of An s ^ fi) {x). One sees that it's 
shape is quite different than any of the other distortions just mentioned as long as one looks over 
a sufficient frequency range, say from 100 — 400 GHz. 

Historically the frequency at which the total S-Z distortion to the spectrum is zero has been 



of great interest. From the eq 11 one may compute the Taylor series describing how this frequency 
changes as the electron temperature increases. Using An^ z 0) (x) and An^ z 0) (x) one may determine 
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Fig. 3. — Plotted versus frequency is ^f- from the classical y-distortion: An^ z 0) (dotted line), 
and from the next order correction: An^ 2 z 0) (solid line). The fact that the solid curve has 

greater amplitude than the dotted curve indicates that ^i 2 ' 0) ~ kT % underestimates the size 

to f Y (l,0) rn-ec 

of this correction. The distinct shape of the correction should allow observations with high 
frequency coverage to separate out this component of the distortion from other sources of foreground 
contamination. 
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the linear term, 

a; (e e ) = 3.83002 + 4.29189^^ + 0(0*) (21) 

which agrees with Fabbri (1981). Rephaeli (1995) has pointed out that this linear extrapolation 
overestimates the shift in xq which suggests that the next term in the series is negative and this is 
indeed found by Challinor & Lasenby (1997). 

Traditionally X-ray measurements are used to determine the gas temperature in galaxy 



clusters. Quantities like those in eq 20 are complementary measures because the X-rays give a 
temperature weighted by the square of the electron density while the S-Z effect is weighted linearly 
with the electron density. Thus X-ray temperatures are more weighted toward the cluster centers 
than S-Z based temperatures. Combining X-ray measurements with S-Z measurements and 
extensions will allow better empirically based modeling of clusters. More generally one can hope 



to make temperature determinations using eq 2C in gas clouds which are too diffuse to produce 



sufficient X-ray luminosity for a good temperature determination. Hubble constant determinations 
must rely on X-ray measurements since here one makes use of the different density dependencies of 
X-ray emission and S-Z distortion. One might consider looking for the finite optical depth effects 
which do depend on a different power of the density, but these are liable to be very difficult to 
measure ( Fabbri 1981| ). 



In summary, it is proposed that one fit for an additional component of the CMBR spectral 
distortion in the direction of galaxy clusters or other hot gas where a y-distortion has been 
detected. The amplitude of this additional component, Y^.o), m ay be used to gain additional 
information about the gas temperature and density. 
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